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We present a cluster dynamical mean-field treatment of the Hubbard model on a square lattice to 
study the evolution of magnetism and quasiparticle properties as the electron filling and interaction 
strength are varied. Our approach for solving the dynamical mean-field equations is an extension 
of Potthoff 's "two-site" method [Phys. Rev. B. 64, 165114 (2001)] where the self-consistent bath 
is represented by a highly restricted set of states. As well as the expected antiferromagnetism close 
to half-filling, we observe distortions of the Fermi surface. The proximity of a van Hove point and 
the incipient antiferromagnetism lead to the evolution from an electron-like Fermi surface away 
from the Mott transition, to a hole-like one near half-filling. Our results also show a gap opening 
anisotropically around the Fermi surface close to the Mott transition (reminiscent of the pseudogap 
phenomenon seen in the cuprate high- Tc superconductors). This leaves Fermi arcs which are closed 
into pockets by lines with very small quasiparticle residue. 

PACS numbers: 71.10.-w, 71.18.-fy, 71.27.-(-a, 75.30.Kz 



I. INTRODUCTION 

Understanding the physics of the Hubbard model con- 
tinues to be a fundamental issue in strongly correlated 
systems. This model captures the transition between a 
metallic state and a correlated insulator; how this tran- 
sition takes place has been investigated by many work- 
ers in the past, with differing approaches emphasizing 
the formation of Hubbard bands, ^ the increasing mass of 
the quasiparticles^ and the proximity to antiferromag- 
netism.- Significant progress has been made in recent 
years by the development of dynamical mean-field the- 
ory (DMFT) 4i Here a set of equations — exact in the limit 
of infinite dimensions — is derived which maps the prob- 
lem onto an interacting impurity model to be solved self- 
consistently. This method has revealed a Mott transition 
that is a synthesis of the pictures of Mott and Brinkman- 
Rice: formation of Hubbard bands in parallel with a 
large mass quasiparticle. The absence of momentum de- 
pendence in correlations within DMFT means that this 
method does not account for variations in quasiparticle 
properties across the Brillouin zone. This is likely to be 
important if the Mott transition has an antiferromagnetic 
aspect. 

Building on the success of dynamical mean-field the- 
ory, extensions to the theory are being actively studied. 
One such extension combines it with density-functional 
theory to improve the treatment of local correlationsi^iSii 
Another considers intersite correlations via so called clus- 
ter DMFTiSi2ii£ It is also often important to include mul- 
tiple local orbitals when modeling real materials such as 
the correlated oxides. In all of these extensions the task 
of solving the resulting self-consistent equations becomes 
increasingly problematic: the equations require the so- 
lution of impurity models with an increasing number 
of local degrees of freedom, which often demand high- 
performance computing resources and sophisticated ap- 



proximation schemes to extract the low energy physics. 

Yet in contrast to these computationally intensive ap- 
proaches, Potthoff^ demonstrated that much of the Mott 
transition physics could be captured with a drastic ap- 
proximation of the non-interacting bath of electrons that 
couples to the impurity in DMFT. Whereas in one com- 
putational scheme for tackling the DMFT equations, the 
bath is modeled by up to twelve coupled sites (the exact 
diagonalization method^^), Potthoff used just a single 
site to represent the bath. The self-consistency condi- 
tions were constructed to ensure that the quasiparticle 
properties and band filling were matched. Solving them 
yielded a successful description of the Mott transition 
showing, for example, the narrowing of the quasiparticle 
resonance and the formation of the Hubbard bands. A 
value of the critical Hubbard interaction, U , was obtained 
comparable to the best calculations. 

In this paper we present an extension of Potthoff's 
"two-site" approach to treat cluster DMFT. We use 
our method to investigate the approach to the Mott 
transition for a single band Hubbard model on a two- 
dimensional square lattice with nearest neighbor hop- 
ping, t. The Hamiltonian is 

k.cr i 

where ek — — 2t[cos(fca;a) -f cos(/cj,a)]. We study this 
model as an example case, although our method is read- 
ily extended to more complicated Hamiltonians. The 
method allows us to quickly investigate the zero temper- 
ature phase diagram across a large range of parameter 
space using a desktop computer. The efficacy of our ap- 
proach is reinforced by our results which are consistent 
with other results in the literature for the magnetic phase 
diagram. Moreover we can go beyond existing work to 
study the quasiparticle properties in momentum space: 
seeing, for example, how an electron-like metal (well away 
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from half-filling) becomes a hole-like doped Mott insula- 
tor near half-filling. Some of our results are suggestive of 
the physics of the cuprate superconductors with the ap- 
pearance of pseudogap regions and "arc-like" Fermi sur- 
faces brought about, in our case, by the combination of 
antiferromagnetism and proximity to a van Hove point. 
We also observe a Fermi surface distortion resulting from 
a Pomeranchuk instability also reported elsewhere in the 
literature ii2ii^4i^*i& Physically, this is due to the proxim- 
ity of the Fermi surface to the van Hove point, but the 
tendency will be exaggerated in our model due to the 
reduced symmetry of our cluster (see later) . 

We begin with a summary of DMFT together with Pot- 
thoff's "two-site" approach. We then discuss how DMFT 
can be extended to a cluster of sites, and describe our ap- 
plication of Potthoff 's approach to a cluster consisting of 
a pair of sites. We then demonstrate this method on the 
2D Hubbard model and present our results. We conclude 
with a discussion of the physics behind these results and 
future extensions of our method. 



II. FROM DMFT TO THE TWO-SITE 
APPROACH 

The DMFT procedure^, can be described as follows. 
We focus on a single site of the Hubbard model, and 
notionally integrate out all the other sites. This gives an 
effective action for the remaining site of the form 

Ses = ~ f dT f dr'cUr)go\r - r')c,(r') 

Jo Jo fr)\ 

.0 ^^-^ 
+ U dTc\{T)cjT)c\{T)c^{T). 

Jo 

The function Qq completely encapsulates the dynamics 
of electrons entering and leaving the site from the rest 
of the lattice; however it is not known a priori since we 
cannot in practice integrate out the other sites. With 
this action we could determine the interacting Green's 
function Giocsdi'i-'^n), and extract a local self-energy from 
the Dyson's equation 

^local(*'^n) = So^i'i-'^n) ^ Slocal(«W„). (3) 

The DMFT ansatz, exact in infinite dimensions, is to 
use this local self-energy as a spatially homogeneous (but 
frequency-dependent) self-energy for the full lattice prob- 
lem: 

Glat {iuJn,ky^ = icJn -f ^ - ek - Slocal (it^n ) • (4) 

The self-consistency requirement that the on-site Green's 
function of the extended lattice (containing the self- 
energy) is the same as the local Green's function we 
started with, i.e. 

Glocal {il^n ) = ^ Glat («W„ , k) , (5) 
k 



provides the constraint on the unknown initial function 
QQ{iujn), thereby completing the self-consistency loop. 

Eq.|21is the effective action of a single interacting impu- 
rity coupled to a continuum bath, but the procedure de- 
scribed above cannot be achieved exactly because the sin- 
gle site impurity problem with an unrestricted bath is still 
intractable. One must approximate and use a model that 
is practically solvable; for example, a finite-sized impurity 
model for exact diagonalizationjiS or a discretized effec- 
tive action for Quantum Monte Carlo methodsii^ This 
means that only limited functional forms of Goii^n) can 
be represented, and also that the final self-consistency 
condition cannot be implemented precisely. 

In the "two-site" realization of DMFT introduced by 
Potthofffii the local model is an impurity site together 
with a bath consisting of a single site only, with a Hamil- 
tonian: 

a 

(6) 

where the electron creation operators a\ and cj. are for 
the impurity site and the bath site respectively. Di- 
agonalizing the non-interacting (C7 = 0) model yields 
Qo(iS) = uj + IjL— [V^ /(w — Cca-)] (where we are now consid- 
ering zero temperature real frequency Green's functions). 
The two-site model allows a minimal frequency depen- 
dence in Qq. By exactly diagonalizing the many particle 
Hamiltonian of Eq. |S1 the local on-impurity-site interact- 
ing Green's function Gi,np(w) can be constructed from 
the Lehmann representation, and the self-energy S(a;) is 
extracted {c.f. Eq. using: 

Grmp,.(^) = ^ + - - Slocal,.(^). (7) 

The full functional self-consistency of Eq. [S] cannot 
be achieved within such a restricted representation so 
it is necessary to decide how best to implement a self- 
consistency requirement. Potthoff chose two physically 
motivated features, taking advantage of the analytic sim- 
plicity of the two-site impurity model. Firstly, the elec- 
tron fillings given by the Green's functions must be 
equal for the impurity model and the lattice model. 
Secondly, features of the central quasiparticle peak are 
matched: the self-energy is reduced to the low energy 
form S]iocai.cr('^) Oct + ba-uj, and terms of the result- 
ing "coherent" impurity and lattice Green's functions 
at high energy are matched (see Ref. for more de- 
tail). In effect, the shape of the central quasiparticle 
peak is analyzed by the size of its high energy tails, iso- 
lated from other parts of the spectrum. The resulting 
self-consistency conditions for the four bath parameters 
(KjEccr) are: 

nimp.cr = niat,o- and = ^ e^Za^ , (8) 

k 
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where the quasiparticle residue 

= ~ -ij$f) ^ (1^ ' 

and we have assumed that X^k — 0. 

Solving these equations produces results for the Mott 
transition which compare well with the full DMFT, and 
properties of the Fermi liquid which are consistent with 
exact results (see Ref. ITTi) . It is hard to imagine a simpler 
model which can do this; it succeeds because of the phys- 
ical motivation of its self-consistency conditions, namely 
the filling and properties of the quasiparticle peak. It 
is a useful approach for calculations on extended models 
such as multiple bands-^ and, as we shall demonstrate, 
clusters. 



III. FROM CLUSTER DMFT TO TWO-SITE 
PAIR-CLUSTER DMFT 

We now consider an extension of this method to clus- 
ter DMFT:^'^ instead of starting with a single site, self- 
consistency conditions are derived for a cluster of sites, 
which allows a momentum dependence in the self-energy. 
The geometry of the lattice becomes important and hence 
different types of magnetic order can be investigated, and 
spectral information varies non-trivially in k-space, un- 
like conventional DMFT. A number of studies have been 
reportedii2ii2i2ai2iiS We describe here a cluster DMFT 
for the case where the cluster consists of just a pair of 
sites. 

We then detail how to implement Potthoff's two-site 
method to solve cluster DMFT. This can be contrasted 
with the exact diagonalization formulation of conven- 
tional DMFTji^ where a single impurity site with multi- 
ple bath sites is used; instead, we use multiple impurity 
sites in a cluster and include correlations in the simplest 
way via Potthoff's scheme for single bath sites. Below, we 
derive the self-consistency conditions for the pair-cluster 
(each of the two cluster sites is connected to a single bath 
site). 

Implementing a cluster DMFT is fundamentally am- 
biguous, as recently noted by Biroli et al.'r^ the approach 



we adopt is "CDMFT" within their classification, which 
is arguably the simplest scheme appropriate for broken 
symmetry states. One of the strengths of conventional 
DMFT is that it is exact in the limit of infinite dimen- 
sions. In contrast, our cluster approach is not a system- 
atic 1/d correction, though it is of course exact in the 
limit of infinite cluster size when the self-consistent bath 
is merely a sophisticated boundary condition.^ 

In cluster DMFT we imagine integrating out all sites 
except those in the cluster. An electron can in general 
now leave from and arrive back at any of the sites within 
the cluster: the function Qq must become a matrix, cou- 
pling together the dynamics of the cluster sites. The 
resulting action (c./. Eq. |2Jl is: 

+ U f dr J2 4,(.rX^ir)4iirKir) , 
ie{A.B} 

(10) 

where the summations are over the cluster sites. Solving 
this local problem now yields a matrix G'iocai,ij and a 
matrix self-energy. 

Different approaches to cluster DMFT involve different 
ways of combining the matrix self-energy with the non- 
interacting lattice Green's function, and different self- 
consistency conditions. We shall now describe our partic- 
ular approach, for the specific case of a 2D square lattice 
with a cluster consisting of a pair of sites. The sites rep- 
resent the two sublattices of the bipartite square lattice, 
and we label them A and B. The impurity model consists 
of a pair of sites connected by a — i hopping element, cut 
out of the original lattice [see Fig. ^a)]. Each site has 
its own independent bath site [Fig.QJb)]. The Hamilto- 
nian for this pair-cluster impurity model is diagonalized 
as previously, resulting in a 2 x 2 matrix Green's func- 
tion, whose self-energy matrix is extracted (c./. Eq. U)) 
using Dyson's equation: 



t-T,AB,a{uj) UJ + fi- i;Bcr(w), 



^imp,.(^) = \ , , y| 1 . J . (11) 



where the matrix index is the cluster site (A, B). 

The non-interacting lattice Green's function now has k-dependence from the Fourier transform of hopping elements 
to surrounding clusters [Fig. QJc)], and is a matrix with respect to the cluster sites {A^B). Combining it with the 
matrix self-energy from the impurity model according to the DMFT ansatz, gives the interacting lattice Green's 
function (as above, c.f. Eq. QJ: 

W('^'^)- l_eke-fc.«-E^s,.(^) c. + M-Es.(^) i' ^^^^ 
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FIG. 1: In our pair-cluster approach we identify a small cluster in the lattice (a), and treat this cluster as an "impurity" 
connected to a bath. A highly truncated basis is chosen for this bath, leading to an impurity model with four sites (b). To 
reconstruct the extended system, the cluster is re-embedded in the lattice periodically (c). The bath is then adjusted to ensure it 
is self-consistent with this embedding. Although our cluster shape manifestly breaks the x-y symmetry, there are self-consistent 
solutions which reflect the original square lattice symmetry. 



Full self-consistency would require that the matrix ver- 
sion of Eq. [SI be satisfied, but here because of the re- 
stricted representation of Go^ij using just two additional 
non-interacting sites, we must adopt a more modest self- 
consistency requirement. We generalize the approach of 
Ref. ^3 and require firstly that the filling of each sub- 
lattice is the same for the impurity and lattice models: 



Mmp,{A,B}.xr — "lat,{A,B},(T 



(13) 



which will effectively constrain ec,{A.B},(7- The second 
condition characterizes the coherent quasiparticle peak. 
Expanding the Green's function matrices at high u, with 
the coherent form of the self-energy in place, gives a ma- 
trix equation self-consistency equation involving a quasi- 
particle residue matrix Z = {I — d'E/dui{0))^^: 



fi t 
t fi 



fi t 
t fi 



E 



-eke 



-Eke" 



-eke 



-eke" 



(14) 



After the k-sum has been carried out we obtain the 
self-consistency condition: 











U^Zac 



(15) 



This constrains {Va,Vb} in terms of the quasiparticle 
residues Z calculated from the local self-energy. 



IV. CALCULATIONS 

We now describe the implementation of our cluster 
DMFT approach. The eight self-consistent equations 



fEasll3landll5|l are solved for the eight bath parameters 
^c,{A,B},ai y{A,B}.cn and additionally we wish to consider 
a given total filling n, and this provides a ninth condi- 
tion which constrains the chemical potential ^ such that 
X] "-{A.Bj.o- = The practical problem is thus nine- 
dimensional root finding, for which we use a Broyden 
method combined with line searches*^ 

The two main computations are diagonalizing the 
Hamiltonian and calculating the band filling for the lat- 
tice case. For the latter, we sum the imaginary part of 
the lattice Green's function (Eq. I12|l over a finite num- 
ber of points in k-space and then integrate numerically 
for oj : (— CX3, 0). A small analytic continuation uj — iS in 
required, with 5 ^ O.OOSt, and typically the k-sum in- 
volves K = 120 points across the Brillouin zone; checks 
were done to ensure that dependence on 5 and K is in- 
significant. Fillings above n = 0.98 are ignored since the 
k-space resolution is insufficient to give an accurate rep- 
resentation of the Fermi surface, except at exactly half- 
filling where this is not an issue. 

For each choice of (J7, n) there is in general more 
than one set of bath parameters which satisfies the self- 
consistency conditions. Phases with paramagnetic, fer- 
romagnetic, antiferromagnetic, ferrimagnetic, insulating 
and charge-ordered character have all emerged. We did 
not consider superconducting order. A pitfall of any self- 
consistency scheme is the production of unphysical ex- 
cited states, and an advantage of the two-site DMFT is 
the ease of calculating the energy of a solution. Using 
Fetter and Walecka^S^ the energy of a fully interacting 
system is 

/dwV- Im tr [fHo(cj,k)-f Js(w)VGiat(t^,k) 

J-OQ ^ Tl^^^-i& W ^ J 

(16) 

where iJo (the non-interacting Hamiltonian) , S and Giat 
are matrices with respect to [A, B) as above. We use this 
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to identify the lowest energy solution. 

As well as the energy and magnetic order, we wish 
to access the k-space spectral information of our solu- 
tions. Combining sublattice Green's functions defined in 



V. RESULTS 

We apply our cluster DMFT approach to the single 
band Hubbard model on a 2D square lattice defined in 
Eq.n We analyze our results for the region < U /t < 30 
and filling 0.75 < n < 1. Wc find a ground state phase 
diagram with a variety of phases, both magnetic and non- 
magnetic, as illustrated in Fig.|21 the extent of magnetism 
is much reduced from that predicted by Hartree-Fock 
theory. Note that all our phases are metallic for n ^ 
1 and insulating for the half-filled case n — 1. From 
an analysis of the full spectral function A{u!,'k) wc can 
characterize the phases and offer a physical origin for 
their stability. 

Recall that for this model there are three important 
factors driving the underlying physics: the Mott transi- 
tion, the nesting of the U — Fermi surface at n = 1, 
and the van Hove points in the free particle dispersion at 
(tt, 0) and (0, tt). Clearly at n = 1 all of these factors play 
a role, but away from this point we see the three factors 
exerting distinct influences on the ground state physics. 



A. Trends across the phase diagram 

The solution far from half-filling (n < 0.82), or for 
very small interaction strength {U < 0.25i), is a metal- 
lic paramagnetic phase [Fig. |21 phase (a)] with a sharply 
defined Fermi surface, seen as a discontinuity in 
[Fig. [TJb), dotted line]. This Fermi liquid becomes in- 
creasingly correlated as U is increased; the quasiparticle 
residue z decreases, and rik becomes more homogeneous 
(though maintaining a discontinuity) as weight spreads 
out through k-space. This reflects the local nature of the 
strong repulsion, U. Spectrally, Hubbard side bands are 
observed to appear, and the central quasiparticle band 
narrows; an example density of states is shown in Fig. O 
consistent with the three peak structure^ used to under- 
stand the Mott transition.A 

As the filling n is increased nearer to 1 (half-filling), the 
Fermi surface in the paramagnetic state expands, and we 
observe a phase transition [to phase (b) of Fig. [2] where 
the Fermi surface breaks x-y symmetry (see Fig. 0J| . It 
is likely that the system is exploiting the reduced Fermi 
velocity near the van Hove point (0, tt) to redistribute 



the reduced Brillouin zone into a single correlation func- 
tion where k is in the extended (non-symmetry-broken) 
Brillouin zone, gives: 



(17) 

I 

electrons in momentum space to lower the total energy — 
a Pomeranchuk transition. This tendency has also been 
seen by other authors li^ji^ii^^i^ Note however, our clus- 
ter breaks x-y symmetry from the outset, so this is not 
strictly symmetry breaking, but the degree of deviation 
from x-y symmetry increases abruptly at this point. 

For U > 6t and n > 0.8, we find that the ground 
state is a spin-symmetry-broken antiferromagnetic phase 
[Fig.|21phase (c)]. The Fermi surface is Pomeranchuk dis- 
torted, remaining the same shape as for the paramagnet 
[Fig. |2] phase (b)]. The local density of states shown 
in Fig. EJa) changes little between that paramagnetic 
and this magnetically ordered state. Indeed, the only 
change is that (n^^ = it-bi) becomes no longer equal to 
{uAi = ns^)- We interpret this phase as a Slater antifer- 
romagnet because the density of states and quasiparticle 
dispersion is, for low U at least, what one would expect 
from a weak spin-density wave transition on a metal- 
lic state with a gap growing like U (see, for example, 
the mean-field theory of Schrieffer, Wen and Zhang^^). 
Note that here the gap in the density of states is already 
present in the paramagnetic state; we will refer back to 
this point in relation to the pseudogap physics of the 
cuprates. 

There is a second antiferromagnetic phase [Fig. [biphase 
(d)] which appears when U > 16t. Here the magnetic 
gap develops within the narrow quasiparticle peak of the 
density of states, which is itself well separated from two 
fully formed Hubbard bands in the familiar three peak 
structure (see Fig. The momentum dependence of 
the narrow quasiparticle dispersion has the same struc- 
tural form as that in the Slater magnetic phase above 
but with a greatly reduced bandwidth. This can be seen 
by comparing the spectral functions shown in Fig. lUb) 
and Fig. El^b) near a; = 0. The magnetism is thus ap- 
pearing like a spin-density wave transition but now for 
the rcnormalizcd quasiparticles. We find that the gap re- 
duces as U is increased consistent with ~ t^/U; thus we 
interpret this phase as an antiferromagnetic metal in an 
effective t-J model, where t^s ^ t and Jcs ^ /U . The 
neighboring paramagnetic phase (at a reduced filling of 
n < 0.9) could perhaps be described, at low energies, as 
a nearly-antiferromagnetic Fermi liquid. 

There are two other magnetic phases. First, exactly 
at half-filling (and extending to lower n near U ~ 5t) 



Gk = 



Lu -\- ^— (Sa(i^) + Ss(w))/2 -I- ek + SyiB(ijj) cos kxa 
{lo + ^i- Syi(w))([j -I- /i - Ss(^^)) - ^^^(cj) - 2Y.AB{uj)e-k cosk^a - 
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FIG. 2: Schematic ground state phase diagram {n,U/t) for 
the square lattice Hubbard model (bandwidth 8t) for our clus- 
ter DMFT, showing a conventional paramagnetic phase (a), 
a paramagnetic phase with a distorted Fermi surface (b), a 
Slater- like antiferromagnetic phase (c), another antiferromag- 
netic phase (d) , a pure Slater nested antiferromagnetic phase 
(e) and a ferromagnetic phase (f). A dashed line indicates a 
second order transition. The phases are characterized further 
in the text. Solutions for the hatched region, 0.98 < n < 1, 
are not considered (see text). 



there is a antiferromagnet [Fig.|21phase (e)] of fairly pure 
Slater character, inevitable for a bipartite lattice. The 
spectrum we find is consistent with that described else- 
where in the literature, for example the separate high 
energy bands in Ref. Second, there is a small patch 
of low moment ferromagnetism [Fig. |21 phase (f)]; here 
the system is exploiting the soft Fermi surface near the 
van Hove points (c./. Ref. Il5j) . and the Fermi surface for 
one spin species distorts more than the other. 
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FIG. 3: At U = 20t, n = 0.8 the ground state is paramag- 
netic [Fig. |5| (a)]. The smoothed local density of states for 
this solution is shown. The central quasiparticle peak shows 
a single-particle density of states, with a bandwidth 8t* ~ 3i 
significantly renormalized from the original 8t due to correla- 
tions. Above and below the peak there are upper and lower 
Hubbard bands, separated by a Mott gap of size ~ U, and 
both Hubbard bands show imprints of the single-particle dis- 
persion. 
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FIG. 4: At i7 = 2t, n = 0.90 the ground state is paramag- 
netic but the Fermi surface has spontaneously distorted [Fig.|21 
phase (b)]. Here the electron filling Wk is shown in grayscale 
across the Brillouin zone: black corresponds to fully occupied 
states, white are completely empty states. The Fermi surface 
is well defined, but has undergone a Pomeranchuk transition: 
electrons have moved from (tt, 0) to (0, tt), exploiting the flat 
dispersion near the van Hove point. Note the appearance of 
a faint "ghost" Fermi surface displaced from the original by 
a wavevector {n /a,n /a). This is due to the effect of antifer- 
romagnetic fluctuations. 



B. Metal— insulator transition and the pseudogap 

A key puzzle in the study of correlated electrons near 
the Mott transition concerns the changes in the Fermi 
surface as the Mott insulating state is approached. How, 
for example, does the large volume Fermi surface ex- 
pected in a Fermi liquid obeying Luttinger's theorem. 
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FIG. 5: At U = 6t, n = 0.88 the ground state is an antiferro- 
magnetic metal with a distorted Fermi surface [Fig. |21 phase 
(c)] . (a) The total local density of states is shown for this state 
(smoothed and combining up and down spins). The magnetic 
gap is apparent, and correlations cause further partitioning. 
In this region the gap size increases with increasing U. The 
densities of states for the neighboring paramagnetic solution 
[Fig- HI phase (b)] are insignificantly different, although the 
sublattice- and spin-resolved densities of states are of course 
unbalanced for antiferromagnetic solutions, (b) The evolution 
of the electron spectral function is shown for various straight- 
line sections through the Brillouin zone. The plot is grayscale 
with darker regions representing states with higher spectral 
weight. As well as weakly dispersing bands at high energies, 
the quasiparticle dispersion can be seen for energies between 
±4i. The form of the dispersion is that expected for a mean 
field spin-density wave state. Hence we identify the gap in 
the density of states as being Slater-like. 



evolve into a doped Mott insulator where it appears 
that the number of holes characterizes physical proper- 
ties (see, for example, Ref. l27j) ? This question can be an- 
swered within the constraints of our approach: as the car- 
rier concentration approaches half-fiUing, we see a large 
Fermi surface breaking up into hole pockets separated 
by regions without free carriers. In this section we will 
describe this process in more detail. 

Well away from half-filling, we observe a distinct Fermi 
surface, defined by a discontinuity in Uk [see dotted line in 
Fig.[7Jb)], with a Luttinger volume equal to n. As ti 1 
and we approach an antiferromagnetic state, evidence of 
a "ghost" Fermi surface is seen: low energy electron-like 
excitations appear at the position of the original Fermi 
surface but displaced by the antiferromagnetic wavevec- 
tor (tt, tt). This can be seen as a faint line in Fig.^ Note 
that the magnetic symmetry has not yet been broken. 

Looking at the complete spectrum [for example in 
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FIG. 6: At U — 20t, n — 0.9 the ground state is also an 
antiferromagnetic metal with little x — y symmetry breaking 
[Fig. |5| phase (d)]. (a) Here we show the local density of 
states (smoothed, up and down spins combined). The central 
quasiparticle peak is separated from upper and lower Hubbard 
bands and is heavily renormalized due to correlations. It is 
also gapped; this gap decreases with increasing U . (b) From 
the spectral density we identify the dispersing quasiparticle 
near the chemical potential. The form of the dispersion is 
very similar to that of Fig.j^b), so Slater-like, but the degree 
of dispersion is very much reduced. 



FigsOc) and (d)] reveals that this "ghost" Fermi surface 
is a manifestation of a "ghost dispersion" — the original 
dispersion shifted by (tt, tt). This ghost dispersion has 
much less spectral weight than the original renormahzed 
quasiparticle dispersion. Nevertheless, it hybridizes with 
the original dispersion to give the final electronic struc- 
ture. The hybridization opens a gap in the dispersion, 
and the effects of the gap depend on whether the chem- 
ical potential /i falls within the gap, and how the gap 
evolves round the Brillouin zone. A schematic view of 
the origin of the electronic structure is given in Fig. |S1 

In the case of Fig. ^ where U — 2t and n — 0.90, 
fi is below the gap everywhere on the Fermi surface, 
and crosses the new dispersion in two places; the spec- 
tral weight at the second crossing is much weaker and 
causes the "ghost" Fermi surface. However, if the filling 
is increased (raising n) or the interaction U is increased, 
widening the hybridization gap, then the chemical po- 
tential may fall within the gap at various places around 
the Brillouin zone. Such an occurrence is illustrated in 
Fig. Eld) for the case {U ^ 5t, n = 0.90) along the Hue 
in the Brillouin zone (0,0) to (tt, tt). With no band cross- 
ing the Fermi surface there is no discontinuity in as 
demonstrated in Fig.jTJb). Instead Uk falls smoothly be- 
cause of the reduced spectral weight in the "ghost" band 
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FIG. 7: At f/ = 5t, n — 0.9 the ground state is paramagnetic and breaks x-y symmetry. It also shows the breakup of the 
Fermi surface into hole pockets by the formation of an anisotropic gap around the Fermi surface. Fig. (a) shows rik in grayscale 
across the Brillouin zone (black <-> 1, white <-> 0). A distinct Fermi surface can be seen along some directions but not others. 
Fig. (b) illustrates the absence of a Fermi surface discontinuity in nk along the line (0,0) —> (tt, tt) (solid line). This should be 
compared to the form of in a phase without such a gap [{U = 8t,n — 0.81) shown by the dotted line]. The nature of the 
anisotropy around the Fermi surface can be seen in the electron spectral functions: (c) A grayscale plot of the spectral function 
(dark larger weight) along the line (0, 0) — > (tt, 0). The chemical potential fi crosses a modified dispersion twice, creating the 
hole pocket — though with little weight at the second crossing point, (d) The spectral function (dark larger weight) along 
the line (0, 0) — > (tt, tt) shows that fi falls within a gap; there is no Fermi surface along this line. The dashed line in (c) and (d) 
shows the non-interacting dispersion. Fig.|H]shows schematically the origin of this electronic structure. 



that forms at the high momentum end of the low energy 
dispersion. There is a gap for single electron excitations 
along this direction. 

This gap is a pseudogap rather than a complete gap 
as there are other directions where a well-defined Fermi- 
surface remains. Fig [7|[c) shows the spectral function 
along the fine (0,0) to (7r,0). There the chemical poten- 
tial crosses the spectral function dispersion in two places. 
The spectral weight is much weaker at the higher momen- 
tum crossing because this part of the dispersion derives 
mainly from the "ghost" band. These two crossings mark 



the edges of a hole pocket as can be seen in Fig. [TJa) . 
The quasiparticles have developed hole-like characteris- 
tics and the Luttinger volume has changed. However, 
because of the evolving spectral weight in momentum 
space, only one side of the pocket has quasiparticles with 
high electron weight. These hole pockets could resemble 
disconnected Fermi surface arcs to experimental probes 
unable to resolve the small electron weight at the second 
dispersion crossing. A similar breakup of the Fermi sur- 
face near the Mott transition was suggested by Furakawa 
et a& 
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FIG. 8: A schematic view of the origin of the structure of 
the electron spectral function in the pseudogap region [as seen 
within our approach in Figs|7[c) and (d)]. Darker lines carry 
more spectral weight, (a) Distinct fluctuations on the A and 
B sublattices result in a second dispersion of low weight (a 
"ghost" dispersion) . This is formed by the folding of the orig- 
inal, possibly renormalized, dispersion in the antiferromag- 
netic Brillouin zone (defining (?m). (b) This hybridizes with 
the original dispersion and forms a gap. Along directions 
in momentum space where the chemical potential lies in the 
gap [line(i)], there is a gap in the excitation spectrum and no 
Fermi surface. Along directions where the chemical potential 
lies away from the gap it crosses the dispersion at two points 
[line (ii)]. A hole pocket forms with low spectral weight at 
the second crossing. 

As n increases, the hole pockets reduce in size until 
they become points, and subsequently vanish and the 
system is an insulator. This illustrates how a smooth 
metal-insulator transition is possible. After a pseudo- 
gap has formed, the Luttinger volume cannot be defined 
as usual; the active Fermi surface consists of hole pock- 
ets, and rest of the filled region of k-space is bounded 
by portions of the original Fermi surface where rik falls 
smoothly, without a step. 

The features described above are reminiscent of the 
cuprate high-temperature superconductors in their nor- 
mal state, which are believed to have a Fermi surface 
consisting of arc-like segmen ts separa ted by pseudogap 
regions (see for example Refs l29l3Ci3lll . A picture where 
the dispersion relation hybridizes with a "ghost" disper- 
sion is consistent with angle-resolved photoemission spec- 
troscopy (ARPES) experiments (see for example Fig. 58 
of Ref. IbOII . and also emerges from other theoretical 
modelsiiS. 

There are important differences, however, between our 
results and what is seen in the cuprates, most signifi- 
cantly in the position of the arcs and the pseudogap re- 
gions. In the cuprates the Fermi surface arcs are seen at 
(7r/2,7r/2) but appear near (0.77r,0) in our calculations, 
and the pseudogap occurs near (tt, 0) in the cuprates but 
near (7r/2,7r/2) in our calculations. However, our results 
are not incompatible, as restriction to a cluster with only 
a pair of sites severely constrains the momentum depen- 



dence of the electron self-energy. Fermi surface shapes 
can only be formed from the original dispersion Ck in 
combination with the asymmetric dispersion cos kxCi 
(as Eg.llTlshows'). This does not admit the d-wave sym- 
metry observed in the cuprates but only s- and p-wave 
configurations, which means for example that arcing can 
only happen near the van Hove points such as (tt, 0). 
Adding more sites to the cluster would progressively opti- 
mize our approximation to the real shape, and we intend 
to extend our work to a four site (2 x 2) cluster which 
would permit d-wave symmetry, though at the price of 
greater numerical complexity. 



VI. CONCLUSION 

We have presented the zero temperature properties of 
a cluster-type extension of dynamical mean-field theory, 
where the self-consistency equations have been solved in 
a very restricted basis. Our approach allows us to calcu- 
late the full phase diagram of a fully interacting system, 
and it can be surveyed and understood through access 
to complete spectral information. The basis we chose is 
"two-site DMFT" , whose results compare favorably with 
other theories and experiment; our cluster extension ad- 
ditionally allows a limited momentum dependence of the 
self-energy. 

Our results provide clues to several puzzling features 
of the metal-insulator transition and how antiferromag- 
netism, proximity to van Hove points, and formation of 
the Mott gap compete with one another; all play a role in 
various regimes. The van Hove point causes Fermi surface 
distortions even at low U; nesting causes Slater antifer- 
romagnetism near half-filling; and the Mott gap becomes 
significant at high U as one might have anticipated. In 
addition, we see how t-J model physics manifests itself at 
large U via the small magnetic gap. We also see how, as 
the metal-insulator transition is approached, the Fermi 
surface evolves from a renormalized Fermi liquid obeying 
Luttinger's theorem, to a pseudogap state where a gap 
opens on some parts of the Fermi surface breaking it up 
into hole pockets (with a strongly momentum-dependent 
spectral density). 

Our work is easily extensible to more sophisticated 
cluster DMFT schemes or multiple bands, and we intend 
to study a 2 X 2 cluster which allows tetragonal symmetry 
and would compare more directly with the cuprates. 
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